Unsharp continuous measurement of a Bose-Einstein condensate: 
full quantum state estimation and the transition to classicality 
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We study a Bose-Einstein condensate (BEC) in a double- well potential subject to an unsharp 
continuous measurement of the atom number in one of the two wells. We investigate the back 
action of the measurement on the quantum dynamics and the viability to monitor the ensuing time 
evolution. For vanishing inter-atomic interactions, mainly the expectation values of the measured 
local observable can be inferred from the measurement record. Conversely, in the presence of 
moderate inter-atomic interactions, the entire many-body state -modified by the measurement- 
is monitored with unit fidelity and, at the same time, the measurement effects a transition from 
quantum to mean- field (classical) behavior of the BEC. We show that this perfect state estimation is 
possible because the inter-atomic interactions enhance the information gained via the measurement. 

PACS numbers: 03.75.Kk, 03.65.Wj, 03.65.Ta, 05.30. Jp 



I. INTRODUCTION 

Continuous measurements enable the monitoring and 
control of individual quantum systems in real time, a 
prerogative in applications of quantum mechanics such 
as quantum information processing and high-precision 
measurements [IHl]. Immanent to quantum physics, the 
trade-off between information gain and disturbance of the 
system requires the continuous measurement to be of low 
eigenvalue resolution per time unit in order to limit the 
measurement-induced alteration of the dynamics. Such 
an unsharp continuous measurement can be realized, e.g., 
by consecutive indirect measurements: the system inter- 
acts in a rapid succession with a sequence of quantum 
probes which are subsequently measured. The strength 
of the continuous measurement can then be controlled via 
the width of the probe's wave function or by tuning the 
interaction with the system. Continuous measurements 
based on sequences of indirect measurements have been 
designed, e.g., to monitor single observables like the po- 
sition of a quantum particle [5 , the photon number [6] in 
cavity QED, the charge distribution in coupled quantum 
dots 7] , and can be used to track Rabi oscillations |8j [9] . 
In further applications, unsharp measurements have been 
employed to estimate the pre-measurement state of an 
ensemble of identically prepared systems [TO] [11] , to de- 
termine the frequency of Rabi oscillations [12], and, com- 
bined with feedback loops, to cool atoms [13] and steer 
the system into a targeted state [14] . 

Recently, ultracold atoms in trapping potentials were 
brought into focus of unsharp measurements as they rep- 
resent interacting many-body systems under exquisite ex- 
perimental control. Among the observations are, e.g., 
the establishment of macroscopic coherence [T5HT7] , even 
in the presence of environmental decoherence [18 , the 
detection and preparation of Mott- or superfluid states 



f!9l [20], the enhancement of macroscopic self-trapping 
[2T] , and the possibility of feedback control [22] . 

In the present contribution, we study a continuously 
measured Bose-Einstein condensate (BEC) in a double- 
well potential. Our motivation is driven by two key- 
targets: On the one hand, we strive to explain the in- 
fluence of the measurement on the quantum many-body 
dynamics. It was shown that the dephasing between the 
BECs in the two wells - which arises due to the inter- 
atomic interactions [23H25] - can be reduced by the mea- 
surement [15] [16] [18] . We unveil the underlying mecha- 
nism and show that the observed behavior is the manifes- 
tation of a transition from quantum to classical (mean- 
field) dynamics of the BEC. 

On the other hand, we address the question whether 
the measurement of a local observable of a many-body 
system can yield complete knowledge of its wave function. 
If so, this would allow complete control of the system 
via feedback depending on the measurement record. To 
tackle this problem, we employ a versatile technique [26] 
which was successfully applied to estimate the state of a 
single particle in several one- and two-dimensional poten- 
tials (among them the classically chaotic Henon-Heiles 
potential) [27] , as well as of a two- level system in the 
presence of noise [28]. We will show, that the strength of 
the inter-atomic interactions in the BEC decides whether 
or not such a state estimation is possible. 

The narrative of this work is as follows: We first intro- 
duce the Bose-Hubbard Hamiltonian in angular momen- 
tum representation as the model for a BEC in a double- 
well potential. In Section [Hi] the realization of the con- 
tinuous measurement of the atom number according to 
Corney and Milburn [16 is discussed together with the 
stochastic equations that describe the coupled dynam- 
ics of state and measurement record. The procedure of 
state estimation is laid out in Section HVl while Sections 
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IVllVIl contains our results: for moderate measurement 
strengths, we analyze quantum dynamics and estimation 
fidelity for the cases of vanishing and weak inter-atomic 
interactions. 



II. MODEL: BEC IN A DOUBLE- WELL 
POTENTIAL 

A Bose-Einstein condensate of N ultra-cold bosons in 
a double-well potential is described by the Bose-Hubbard 
(BH) Hamiltonian [23l [29] : 

H BH = UJ 2 Z -KJ X , (1) 

with angular momentum operators J x — (&{&2 + bib\)/^, 
J y = ifyh - b\b 2 )/2, J z = (b\h - b\b 2 )/2, and J+ = 
b\b2. Here, b^ are the bosonic annihilation (creation) 
operators, and hi = b\bi is the number counting opera- 
tor at site i [30] . In ([!]), U and K parametrize the on-site 
inter-atomic interaction and the tunnelling strength, re- 
spectively. Experimentally, both parameters can be inde- 
pendently controlled via the height of the potential bar- 
rier and by additional magnetic fields that induce Fesh- 
bach resonances [31j[32]. Apart from the total energy E, 
also the total particle number n\ -\-ri2 = N = 2j is a con- 
stant of motion since [J 2 , Hbh] = 0. Thus, the states of 
the Fock number basis, i.e., the eigenstates of J Z: can be 
equally expressed as 1 711,712) = \ n i, N — ni) = |m), where 
the angular momentum m = n\ — N/2 ranges from —j 
to j = N/2. The physical interpretation of the opera- 
tors Ji is as follows: J z measures the particle imbalance 
between the wells, while J y represents the condensate's 
momentum, and J x bears direct information about the 
relative phase of the condensate's fractions in the left and 
right well. The corresponding expectation values (•) with 
respect to the quantum state are readily expressed with 
the help of the single-particle Bloch- vector 
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S= (S X , S y , S z ) = — ((J x ), (Jy), (Jz))- (2) 

In the mean-field or classical limit of large particle 
numbers N (at fixed value UN) the dynamics of the 
condensate is described by the discrete Gross-Pitaevskii 
(GP) equation. The main assumption underlying the 
latter is that the quantum state remains, at all times, 
a SU(2)- or atomic coherent state, i.e., a state of mini- 
mal uncertainties, with respect to the components of an- 
gular momentum (see, e.g., [23j[24]). The coherence of 
the evolving state is sensitively measured by the single- 
particle- (or one-body) purity 

P = ^(1 + 4 + 4 + ( 3 ) 

which ranges from 1/2 to 1 and takes the maximal value 
of one only for coherent states. Thus, a value of p < 1 
indicates the departure from the mean-field, or classical, 




FIG. 1. (color online). Continuous unsharp detection of the 
atom number 712 in the right well of a double-well potential: 
the right well of the atom trap (red) is placed in a leaky cav- 
ity, which is illuminated by a continuous-wave laser detuned 
from atomic transitions. Due to the atom-light scattering 
proportional to 712, a relative phase shift is imprinted on the 
laser light which is subsequently detected by a Mach-Zehnder 
interferometer and read off as the photo current i(t). 

behavior. Further information on the state is contained 
in its Wigner function pw, a quasi-probability distribu- 
tion evaluated on the spin-j Bloch sphere as a function of 
the polar and azimuthal angles and (/>, respectively (see, 
e.g., [33] and references therein). For coherent states, the 
pw take a Gaussian shape, in particular they are posi- 
tive. A coherent state can thus be described by its center 
alone (a c- number), the evolution of which is given by the 
GP equation. In angular momentum representation, the 
latter can be expressed with the help of the Bloch vector 

s x u SyS z1 Sy = s z -\- u s x s zi s z = Sy : (4) 

where the time in Q is rescaled by t —> tK and de- 
scribes the mean-field dynamics on a spin-j Bloch sphere 
which is governed by the control parameter u = UN /K 
[23| l34| 135] . Throughout the paper, we focus on the 
so-called Rabi regime of weak inter-atomic interactions 
u < 1, where the atomic collisions do not yet induce 
self-trapping [36] but notably influence the quantum dy- 
namics [23-25 of the condensate, as the latter oscillates 
between the two wells. Quantum and mean-field dynam- 
ics of the un-monitored system are revised further down 
in Sections IVl and [VTl 



III. CONTINUOUS MEASUREMENTS 

Several proposals outlined the non-destructive measur- 
ing of a BEC via light- [TSHHJ [37] or particle scattering 
[38H4Q] including feedback control [22 . Here, we focus 
on the scheme of Corney and Milburn [16] sketched in 
Fig. [l] In this setup, the number of atoms in one, say 
the right, site of a double- well potential is continuously 
and unsharply measured by placing the corresponding 
site inside a leaky optical cavity. Well-detuned light from 
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a continuous-wave laser pumps the cavity and is not ab- 
sorbed by the atoms but suffers a phase shift proportional 
to the number of present atoms, which is detected 
by superposing the output light with a reference beam 
in a Mach-Zehnder interferometer (homodyning) . Also 
the atoms in the well experience a phase shift induced by 
the AC- Stark effect which is proportional to (J z ) and can 
thus be compensated by tilting the double-well potential 
in the earth gravitational field [IB] . 

After adiabatic elimination of the degrees of freedom 
of the light field, the master equation for the reduced 
density operator p of the BEC reads [T6j [41] : 



-iU[jlp]+iK[J x ,p} + J[J Z ,[J; 



(5) 



with h = 1. The strength of the measurement 7 = 
64% 2 e 2 /r 3 is expressed in terms of the average inter- 
action strength \ between atoms and light field in the 
cavity, the amplitude e of the coherent field pumping 
the cavity and the cavity damping rate T. For the mea- 
surement strength of 7 = 0.01K used below, the optical 
field strength amounts to 0.75mW, where the mass of the 
bosons was assumed to be 1.5 x 10 _25 kg [IB] . 

The master equation (|5| describes the evolution of the 
mixed state associated with an ensemble of BECs which 
are prepared initially in the same state but then undergo 
evolutions with different measurement records, i.e., dif- 
ferent detected photo currents (cp. setup in Fig. [I]). The 
first two terms in ([5| correspond to the unitary evolu- 
tion given by the Bose-Hubbard Hamiltonian ([I]) while 
the third term represents decoherence due to the interac- 
tion with the light beam. This decoherence occurs with 
respect to the eigenbasis of the measured observable. As 
the total number N of atoms is fixed in this setup, a mea- 
surement of the number of atoms in one well is equivalent 
to a measurement of J z . Thus, the third term is the sig- 
nature of the continuous measurement of the number of 
atoms in one well. 



A. Selective regime and stochastic Schrodinger 
equation 

In the following, we do not consider the mixed state 
that would result from averaging the possible outcomes 
of the continuous measurement. Instead, we are inter- 
ested in the state evolution given a particular measure- 
ment record, i.e., a particular time-dependent photo cur- 
rent i(t) induced in the photodetector (cp. Fig{I]). In this 
so-called selective or conditional regime of measurement, 
i.e., conditioning the state evolution on a given photo 
current i(t), an initially pure state \ip(to)} of the conden- 
sate stays pure and can be described by a (stochastic) 
Schrodinger equation [IB] : 



d\M*)) 



[~iHBH-^(j z -{J z )c) 



dt\Mt)) 



where Hbh is the Bose-Hubbard Hamiltonian |l]) and 
(J z ) c is the expectation value with respect to the state 
\ip c (t)). The subscript c marks expectation values and 
system states that are conditioned on a measurement 
record, i.e, conditioned on a specific photo current i(t) 
measured until the time t. Instead of i(t) (cp. Fig{l]), 
we consider the (proportional) signal I(t) := 2i(t)/j : to 
keep our notation simple. Then the increment of the 
measurement signal I(t) from time t to t + dt reads: 



dl= (J z ) c dt + j- 1/2 dW , 



(7) 



where dW is the increment of a white Gaussian noise pro- 
cess (so-called Wiener increment) [42 , that reflects the 
measurement noise. As in the master equation ([5|, the 
first term in the stochastic evolution (|6| represents the 
Hamiltonian (unitary) dynamics while the second (non- 
linear) term, corresponds to the double commutator in 



5|. As it is diagonal in J z , the latter term does not al- 
ter (J z ), but reduces the expectation values of J x and 
J y . The third, stochastic term -in combination with the 
second one- narrows on average the variance of the mea- 
sured observable J z . To illustrate the action of the un- 
sharp measurement with a simple example, we remark 
that in the absence of Hamiltonian dynamics, a continu- 
ous measurement projects the system (due to the second 
and third term of (6|) asymptotically into an eigenstate 
of the measured observable and this eigenstate can be 
identified by the measurement record [43] . 

In the present case of a BEC evolving in a double-well 
potential, different energy scales compete. From both, 
the master equation ([5| and the stochastic Schrodinger 
equation (J6|, it is evident that terms corresponding 
to measurement and inter-atomic interactions increase 
quadratically with N while the tunneling coupling de- 
pends linearly on N. It is thus sensible to consider 
7 := ^N/K as the normalized measurement strength 
with respect to the inter-well tunneling (accordingly, U 
should be directly compared to 7). The larger the value 
of 7, the faster individual occupation levels ri2 can be 
resolved due to the measurement [8, 44 . For 7 ^> 1, this 
results in the so-called shelving, i.e., while tunneling, the 
system remains longer on the individual levels 77.2 an d the 
tunneling oscillations are distorted [45]. In our present 
contribution, we concentrate on small to moderate mea- 
surement strengths 7. The exemplary value of 7 = 1, 
used later, fulfills the operative definition that the tun- 
neling frequency of the non-interacting system (u = 0) 
be barely changed, i.e., that we are far away from the 
shelving regime (see Sec.VB further down). 



;(Jz ~ (Jz)c)(dl - (j z ) c dt)\Mt)) > (6) 



As a side remark, let us note that similar physics, as de- 
scribed by the above master equation ([5|, arises when the 
noise in the quantum dynamics does not result from an 
unsharp measurement but from other sources like, e.g., 
interactions of the BEC with non-condensed atoms or by 
enforced stochastic driving 
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IV. STATE ESTIMATION 

According to Eq. (|7|), the measurement record i(t) con- 
tains information about the expectation value of the mea- 
sured observable J Z1 which can be extracted for suffi- 
ciently high measurement strength 7 using standard tech- 
niques such as Wiener filters [50], see, e.g., [8 . Here, we 
would like to proceed beyond this point, namely, to infer 
(estimate) from a local observable J z the entire quan- 
tum many-body state as it dynamically evolves under 
the influence of the measurement. Moreover, this is to 
be accomplished in a single run, not in a state tomog- 
raphy experiment [5lJ [52] which requires repeated mea- 
surements on equally prepared systems. We will assume 
that only the Bose-Hubbard Hamiltonian is known (i.e., 
the total particle number N and the control parameter 
u) but not the initial state |V^o)) of the condensate. We 
then continuously update an initial guess of the state 
|^e(to)) according to the measurement results (i.e., the 
photo current i(t)) in order to obtain an estimate \ip e (i)) 
of the real state \ip c (t)). The overlap between both states 
is a measure for the fidelity of the estimate: 

F(t) = \(Mt)\Mt))\ 2 ■ (8) 

For the special case of a continuously measured par- 
ticle that maintains a Gaussian-shaped wave function, 
state estimation was discussed in [53j [54] . The estima- 
tion scheme we consider here [26] is versatile and can be 
employed in any continuous measurement scheme. It is 
based on the Ito-formalism and for ideal [26] continuous 
measurements of otherwise closed, quantum systems, an- 
alytic arguments have been given that indicate the con- 
vergence of the estimate to the real state except for cer- 
tain marginal cases [26] (for a recent discussion on this 
convergence see, e.g., [55, 56 ). In this scheme, the state 
estimate is propagated in the same way as any real state 
conditioned on a given measurement record I(t): 

d\Mt)) = (-iHbh ~\{Jz- (Jz)e) 2 ^j dt\Mt)) 

+ \{Jz ~ (Jz)e)(dl ~ (J z ) e dt)\Mt))- (9) 

At first sight, this equation appears to be identical to (|6| 
with the indices c and e exchanged. Note, however, that 
the measurement signal I(t) varies about the expectation 
value of J z with respect to the real state rather than to 
the state estimate. This reflects the fact that \ip e (t)} 
is updated according to the measurement current I(t) 
which is obtained with the system being in the real state. 
Therefore, the state estimate is slaved (i.e., linked) via the 
measurement signal I(t) to the real state evolution. 

V. RESULTS I: VANISHING INTER- ATOMIC 
INTERACTIONS 

We now have the tools at hand to study the effect of 
the unsharp measurement on the system dynamics and 



vanishing interactions u=0 




FIG. 2. (color online) System dynamics for vanishing inter- 
actions u = Upper panels: a) The black lines indicate the 
mean-field phase space of a BEC in a double-well potential 
given by the discrete Gross-Pitaevskii equation Q for vanish- 
ing interactions u = in the absence of measurement 7 = 0. 
The red line marks the evolution of the quantum mechanical 
Bloch vector s for N = 100 particles, located initially in the 
left well (North pole of the Bloch sphere) and times t = — 10 
Rabi periods The Bloch vector remains (for all times) on 
the surface of the Bloch sphere, i.e., the state remains coher- 
ent throughout the evolution, b) The corresponding Wigner 
function pw after w 4 Rabi periods which remains Gaus- 
sian and positive (red areas) [57 . Middle panel: The black 
(blue) line denotes the population imbalance (J z ) c ({J z )e) of 
the real (estimated) wave function. The red line indicates the 
fidelity F(t) while the green line corresponds to the one- 
body purity pjt) . All quantities are plotted in multiples of the 
Rabi period = 2tv/K. Lower panels, d) Same as panel 
a) but now for non- vanishing measurement strength 7=1 
and times t — 10 — 40 tR. The Bloch vector spirals towards 
the center of the sphere and is smallest at about t = 33£i? 
where also the one-body purity (see panel c) assumes its min- 
imal value, e) The corresponding Wigner function pw after 
« 25 Rabi periods (marked by a black arrow in panel c): A 
strong non-coherent structure due to quantum interferences 
has emerged, in stark contrast to the case 7 = shown in 
panel b). 
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FIG. 3. (color online) Details of the system dynamics for 
vanishing interactions u — Upper panels: As in Fig. ^p), the 
black solid (blue dashed) lines denote the angular momentum 
expectation values (-) c ((-) e ) with respect to the real (esti- 
mated) wave function. The time is plotted in multiples of 
the Rabi period tR = 2tt/K. While (J y ) e and ( J z ) e approach 
the values of the corresponding real wave function \ijj c (t)) : the 
estimate (J x ) e evolves stochastically. Lower panels: The vari- 
ances corresponding to the quantities in the upper panels are 
shown. 



to explore whether we can successfully estimate the state 
of the entire system from the record of the local measure- 
ment on J z alone. At the focus of our work is how an 
initially coherent state is affected by the measurement, 
rather than the measurement- induced buildup of coher- 
ence studied, e.g., in [T5HT7] . Hence, we concentrate our 
efforts on the exemplary case of all bosons being initially 
located in the first well, i.e., \i/j c (t = 0)) = an exper- 
imentally routinely prepared SU(2)-coherent state [58] . 
located at the North pole of the Bloch sphere [59] . 

As for the initial state of the estimate \ip e (t = 0)) = 



E 



|m), we pick the "maximally uncertain" state, 



i.e., an equally weighted superposition of 0,1,..., iV 
atoms in the first well |c m | 2 = [2j + l] _1 with randomized 
relative phases arg(c m ). Throughout all our numerical 
calculations, we fix the boson number TV = 100 and the 
inter- well tunneling coupling K = 1. 



weak interactions u=1 
a) b) 






FIG. 4. (color online) System dynamics for weak interactions 
u — 1 Upper panels: a) The black lines indicate the mean- 
field phase space of a BEC in a double-well potential given 
by the discrete Gross-Pitaevskii equation Q for weak inter- 
actions u = 1 in the absence of measurement 7 = 0. The red 
line marks the evolution of the quantum mechanical Bloch 
vector s for N = 100 particles, located initially in the left 
well (North pole of the Bloch sphere) and times t = — 10 
Rabi periods. The Bloch vector spirals into the interior of the 
Bloch sphere, indicating a loss of the one-body purity, b) The 
corresponding Wigner function pw after t = 10t# develops 
an intricate structure with strong positive (red) and negative 
(blue) contributions [57 . It is also smeared out along the cen- 
tral classical orbit of panel a) which results in s y — s z — 0. 
Middle panel: The black (blue) line denotes the population 
imbalance (J z ) c ({Jz)e) of the real (estimated) wave function. 
The red line indicates the fidelity F(t) (I8J) while the green 
line corresponds to the one-body purity p(t). Lower panels, 
d) Same as panel a) but now for no n- vanishing measurement 
strength 7=1 and times t = 30 — 60 tR. The Bloch vector es- 
sentially remains on the surface of the sphere and follows the 
mean-field (classical) evolution. The corresponding one-body 
purity stays approximately at one, in stark contrast to the 
case 7 = shown in panel a), e) The corresponding Wigner 
function pw after « 48 Rabi periods (marked by a black arrow 
in panel c): The latter remains essentially a coherent state, 
opposite to the case 7 = (panel b). 
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A. Un-monitored dynamics 

We start our discussion with the special and instruc- 
tive case of vanishing inter-atomic interactions u = 
which was barely touched upon in [16] and for which we 
first revise the un-monitored dynamics. Without mea- 
surement 7 = 0, the Bloch vector s of the system simply 
rotates at the Rabi frequency K on the meridian defined 
by the intersection of the Bloch sphere with the y — z 
plane [24], as shown in Fig. . In this Rabi regime, 
the Wigner function pw associated with the quantum 
state, retains its Gaussian shape (cp. Fig. J^)) and its 
center moves on the corresponding mean-field trajectory, 
i.e., the respective solution of the GP equation Q for 
u = [57]. Together with the lack of negative values, 
this points out a "classical state" which possesses a joint 
probability distribution of the angles and <ft and, hence, 
state estimation might be expected to work best for this 
case. 



B. Monitored dynamics 

The effect of the measurement on the quantum dy- 
namics is shown in Figs, [^-d) for the case 7 = 1. From 
panel c), we read off that the oscillation period of the 
population imbalance (J z ) c remains approximately the 
Rabi-period tn = 2tt/K (in accordance with our opera- 
tive definition of a moderate measurement strength, see 
Sec. III). It is also evident, however, that the amplitude 
of the oscillations is not constant (as it is for 7 = 0) but 



fluctuates. The decrease in oscillation amplitude is ac- 
companied by a reduction of the one-body purity p(t) (cp. 
Fig. [2^)) which indicates that the wave function looses its 
Gaussian shape, i.e., departs from the mean-field behav- 
ior. This is further corroborated by the Wigner function 
of the real state \ijj c (t)) taken at about t « 25£# (cp. 
Fig. [2^)) in which a prominent structure of both posi- 
tive and negative contributions is evident. Furthermore, 
while pw stays well localized with respect to the polar 
angle (i.e., a small variance in there is no permanent 
localization in the perpendicular, i.e., in the azimuthal in- 
direction. This spread over a part of the corresponding 
great circle on the Bloch sphere leads to a reduction of the 
length of the Bloch vector s c and, thus, degrades p(t) and 
the amplitude of the oscillations of (J z ) c - Summarizing 
these observations, one may thus say, that for vanishing 
inter- atomic interactions u = 0, the measurement steers 
the system away from the mean-field behavior revised in 
the previous paragraph. 

We present two complementary perspectives on this 
phenomenon, which we also observed for other ratios 
7 < 1 of measurement and inter- well tunneling strengths: 
in the Schrodinger picture, the unsharp measurement of 
the population imbalance reduces the variance of the cor- 
responding observable J z (cp. Fig. |i]f)) due to the gain 
of information on its expectation value. As a result of 
the unitary dynamics (which amounts to a bare rota- 



tion about the x— axis for u = 0) this also yields a (sys- 
tematic) reduction of the variance of J y (cp. Fig. K3e)) 
but not of J X1 i.e, it does not give rise to localiza- 
tion in the x— direction (cp. Fig. |3]i). Seen in the 
Heisenberg picture, the unsharp measurement reduces 
the variance of the time-dependent observable J z (t) mea- 
sured over a time period At. In the case of vanishing 
atomic collisions u = 0, the orbit of J z (t) is given by 
{s'm(Kt)J y + cos(Kt)J z \t G At}, due to the rotation 
about the x— axis. It does not include the angular mo- 
mentum component J x (the latter is constant and, in our 
case, its expectation value equals zero during the time 
evolution). Thus, the variance in J x is not (systemati- 
cally) reduced but evolves in a predominantly stochastic 
fashion. In the above example, the variance in J x in- 
creases up to t « 33^, where p(t) assumes its minimum, 
and then decreases again (cp. solid line in Fig. [3]i)). 

As a somewhat crude mechanical analogy of this pic- 
ture imagine a potter's wheel that spins about the x— axis 
and the chime representing the uncertainty of the quan- 
tum state in the three directions. The effect of the un- 
sharp measurement is like softly compressing the chime 
in z— direction: due to the wheel's rotation, one will fi- 
nally obtain a piece of chime that is cylindrical, with 
small radius (i.e., uncertainty) in the y — z-plane but 
increased height in the x— direction. One aspect which 
is not captured by this mechanical picture is that, due 
to the stochastic character of the unsharp measurement, 
there is an additional random component in the evolu- 
tion of the chime's shape. That is, the uncertainty in J x 
is not monotonically growing and the one in J y and J z is 
not monotonically decreasing. 



C. State estimation 

We have yet to discuss the performance of the esti- 
mator ([9| to infer the many-body state of the system 
from the measured photo current z(t). From the exem- 
plary case shown in Fig. [2j we find that the oscillations in 
(J z ) c themselves are, after a short transition time, well 
reproduced by the estimator (J z ) e . Nevertheless, the es- 
timated state \^ e (f)) does not fully converge to the real 
state \ip c (t)) within the considered time period of 60 Rabi 
periods, as spelled out by the fidelity that stays below 
unity (cp Fig. [2^)). Given the fact that the un-monitored 
BEC dynamics is "simple" Rabi oscillations in the popu- 
lation imbalance, this is somewhat surprising. Moreover, 
the fidelity oscillates quite strongly, despite the rather 
good agreement of the real (J z ) c and the estimated (J z ) e 
particle imbalance. An inspection of the estimated state 
reveals (cp. Figs. [3^l)-c) ) that indeed (J y ) e an d (J z )e 
agree quite well with the corresponding components of 
the real state while (J x ) e behaves erratically. We can 
immediately understand this behavior by means of the 
above discussed evolution in the Heisenberg picture: due 
to the rotation about the x— axis, the information on J x 
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is not available in the measurement and, hence, the esti- 
mated (J x )e does not converge to the real result (J x ) c . 



VI. RESULTS II: WEAK INTER- ATOMIC 
INTERACTIONS 

A. Un-monitored dynamics 

We now turn to the more general scenario of non- 
vanishing but weak inter- atomic interactions u < 1. In 
this case, the mean-field (i.e., the classical) approxima- 
tion deviates substantially from the quantum dynam- 
ics, already in the absence of measurements (see, e.g. 
[23] [24]): the mean-field solutions are again closed trajec- 
tories on the Bloch sphere with pronounced oscillations 
in the particle imbalance s z (see black lines in Fig. |4^l)). 
Due to the inter-atomic interactions, the dynamics is not 
merely a tunneling-induced Rabi-like rotation about a 
fixed (here the x~) axis but there is an additional "non- 
linear rotation" about the z axis with a rotation angle de- 
pending on the ^-component (see Hamiltonian 0) and 
a different rotation frequency for each trajectory [60] . 
Since the quantum state for a finite number N of atoms 
is not a single point on the Bloch sphere, but has a 
width according to its Wigner function, the different evo- 
lution frequencies lead to a dephasing in the course of the 
tunneling dynamics. As a result, the oscillations in the 
atomic population imbalance are damped and cease af- 
ter a few cycles, even in the regime of weak inter-atomic 
interactions u = 1. In this process, the quantum mechan- 
ical Bloch vector s spirals onto the x— axis of the Bloch 
sphere which implies that eventually (J y ) = (J z ) = 
(see red line in Figjik)) [61]. The corresponding Wigner 
function (cp. Fig. Hp)) assumes positive (red) and nega- 
tive (blue) values and, after several Rabi periods is 
almost symmetrically distributed about the x-axis which 
again implies (J z ) = 0, and thus the collapse of the os- 
cillations. 



B. Monitored dynamics 

It may come as a surprise that the impact of the un- 
sharp measurement of J z is to restore the oscillations, 
even for moderate strengths 7 = 1 (cp. Fig. [4j:)) [62] . 
This effect was first mentioned in [I6j [18] and we would 
like to add to its understanding with a detailed analysis 
supplemented by a mean-field perspective before we pro- 
ceed to the entirely new results on the state estimation. 

We start with the observation that in the presence of 
inter- atomic interactions u ^ 0, the unsharp measure- 
ment induces classical behavior of the BEC: (i) in con- 
trast to the un-monitored quantum dynamics for u = 1 
discussed in the beginning of this section, under the un- 
sharp measurement (e.g., 7 = u = 1), the Bloch vector 
s c essentially stays at the surface of the Bloch sphere 
(cp. Fig. ) with p(t) w 1 (cp. green line in Fig. [4^)). 



That is, the BEC stays approximately in a coherent state 
of minimal uncertainty with respect to the angular mo- 
mentum components. This is further corroborated by 
our observation, that the Wigner function of the BEC 
is positive and approximately Gaussian (cp. Fig. 4^)). 
To draw the connection to the above discussed mechani- 
cal analog, the spinning axis of the potter's wheel is now 
changing in time, due to the nonlinear rotation caused by 
the nonlinear term in ([I]). As a result, the variances of 
all three angular momentum components are systemati- 
cally reduced (still subject, however, to fluctuations that 
result from the stochastic character of the measurement). 

This has to be contrasted with the case of vanishing 
interactions u = 0: here the un-monitored dynamics fol- 
lows the classical evolution with p(t) = 1 but once the 
measurement is turned on, becomes truly quantum as 
indicated by a non-positive Wigner function and a fluc- 
tuating one-body purity p(t). 

(ii) Secondly, we find that this approximately coherent 
state evolves in the vicinity of the classical (mean-field) 
trajectories (cp. Fig. ) . In other words, by means of 
continuous unsharp measurements, classical trajectories 
can be realized, that would rapidly dephase for 7 = 0. In 
this sense, the unsharp measurement can be considered 
weak -as it does not considerably change the mean-field 
dynamics. But, at the same time, it has a pronounced 
impact on the quantum dynamics, namely, to halt the 
dephasing. 

Let us, however, point out that, as a result of the 
stochasticity introduced by the measurement, the quan- 
tum state does not evolve on a single mean-field trajec- 
tory, i.e., a line on the Bloch sphere, but explores an 
entire area (cp. Fig. [4]i)). As a side remark, we note 
that this can as well be inferred from the time evolu- 
tion of (J z ) c shown in Fig. [i^): as mentioned above, in 
the presence of interatomic interactions u 7^ 0, differ- 
ent mean-field trajectories come with different frequen- 
cies [33]. We should thus observe that the oscillation in 
the particle imbalance (J z ) c are restored [16] [18] but not 
necessarily at the same period. Indeed, we find that the 
frequency of the oscillations in Fig. [4^) changes during 
the course of time (see, e.g., around t = 50^), while this 
is not the case for u — (cp. Fig. (2^)). We stress, that 
the smaller oscillation amplitude around t = 50tR corre- 
sponds to an orbit located at the s x < hemisphere, as 
can be seen from Fig. [4]i). Specifically, it is not due to a 
dephasing as for u = since, in contrast to Fig. [2^), the 
one-body-purity stays close to one. 



C. State estimation 

After discussing the underlying mechanism responsi- 
ble for halting the dephasing, we now turn to the perfor- 
mance of the estimator In contrast to the case u = 0, 
not only are the oscillations of (J z ) c well described by es- 
timate (J z ) e but also the fidelity increases within a few 
cycles to unity (cp. Fig. [4^)), indicating the convergence 
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of the estimated state \ip e (t)} to the real state \ip c (t)). 
That is, a local measurement of the particle imbalance J z 
leads to information on the entire wave function \ip c (t)}. 
To understand this, we return to the Heisenberg picture, 
discussed above in Section [V] as a result of the nonlinear 
rotation for u ^ 0, the set of measured observables is no 
longer restricted to the y — z-plane but also contains J x . 
Hence, the measurement of J z yields direct information 
on all angular momentum components and, thus, we ob- 
tain an informationally complete [63] set of observables. 
That is, sufficient information can be gathered from the 
continuous unsharp measurement to determine the state 
of the system. Indeed, without atomic collisions (u = 0) 
where the orbit of J z (t) does not contain J X1 the esti- 
mate does not converge towards the real state within the 
considered time period. 

We emphasize, that the complete knowledge of the 
wave function was gained via the measurement of a single 
observable on a single realization of the system consist- 
ing of hundred atoms with a Hilbert space of dimension 
d = 101. This is remarkable because a full state tomog- 
raphy (by means of von Neumann projection measure- 
ments) of such a system would require measurements of 
d 2 — 1 = 10200 linear independent observables, each one 
carried out on a large ensemble of equally prepared sys- 
tems |64]. 

VII. SUMMARY AND DISCUSSION 

In summary, we studied the continuous unsharp mea- 
surement of the atom number in one site of a double- well 
potential loaded with a BEC. Based on the setup pro- 
posed in [16 we focused on the selective regime of mea- 
surement. In a first step, we investigated the interplay 
between internal many-body dynamics and the measure- 
ment for the case of an initially coherent state of maxi- 
mal particle imbalance. We showed that in the absence 
of inter-atomic interactions, an unsharp measurement of 
moderate strength steers the system away from the coher- 
ent mean- field evolution of the un- monitored BEC, and 
leads to a genuine quantum state manifest in a Wigner 
function of oscillating sign. Contrary, for weak inter- 
atomic interactions (u = 1) the un- monitored dynamics 
rapidly deviates from the mean-field behavior (23[ [24] 
and the inter-well tunneling oscillations of the BEC are 
damped out. As first noted in p~6j [18], the latter can 
be restored by a measurement of moderate strength 7. 
We explained the underlying mechanism and showed that 
the monitored dynamics evolves essentially as a coherent 
state, i.e., close to the surface of the Bloch sphere and 
in the vicinity of the solutions that correspond to the 
Gross-Pitaevskii equation. Hence, we demonstrated that 
the unsharp measurement can induce quantum or clas- 
sical (mean-field) behavior depending on the strength of 
the inter- atomic interactions. 

Secondly, we went beyond the standard measurement 
scheme and explored the viability to infer from the mea- 



surement record not only the expectation value of the un- 
sharply measured observable but the entire system state 
[26] in a single run of the experiment and without any 
knowledge of the initial state. Opposite to the naive ex- 
pectation, this was not achieved for u = but in the 
presence of weak inter-atomic interactions u 7^ and suf- 
ficiently high but still moderate measurement strength 
7. Specifically, we demonstrated for the exemplary case 
of u = 1 that, after a certain waiting period, the time- 
evolving state of the atoms can be monitored with perfect 
fidelity. This convergence in the course of a continuous 
measurement was explained in terms of an information- 
ally complete set of observables that arises from the in- 
teractions. 

The observations for u 7^ 0, show that the continu- 
ous measurement of moderate strength 7 = 1, modifies 
the dynamics considerably as it compensates the dephas- 
ing and induces classical (mean-field) dynamics. On the 
other hand, it does not alter the structure of the mean- 
field dynamics itself. That is, combined with our estima- 
tion procedure, the unsharp measurement can be used 
to prepare a non-dispersing coherent wave packet that 
evolves close to the mean-field solutions and, at the same 
time, monitor the state with perfect fidelity. This should 
be contrasted with recent studies on periodically [48] or 
stochastically [46] [49] driven BECs in double- well poten- 
tials, where due to a strong driving, the dephasing of an 
initially coherent state is reduced but at the same time, 
the entire dynamics is brought to a halt. It is as well dif- 
ferent from the dynamical localization observed for atoms 
in periodically driven lattices (see, e.g., [65]) where by 
appropriate choice of the driving parameter the inter-site 
tunneling is effectively switched off. As a last example for 
the mitigation of dephasing, we mention non-dispersing 
wave packets that have been observed in the context of 
Rydberg atoms [66H68] and which may actually exhibit 
dynamics. In that case, however, the dephasing is sup- 
pressed due to an external driving which typically alters 
the underlying classical phase space considerably. 

Let us discuss the convergence of the estimate and the 
emergence of classical properties from the perspective of 
measurement theory. The arguments given in [26 , sup- 
ported by our numerical evidence for u ^ 0, indicate that 
the convergence of the estimated state to the real state 
-conditioned on a particular measurement record I(t)- 
occurs independently of the choice of the initial estimate 
\ip e (t — 0))- If so, also all real initial states conditioned 
on the same I(t) converge to the same state \^ c (t)) given 
by the measurement record since estimated and condi- 
tioned states propagate according to the same evolution 
equation, cp. Eqs. ([6| and Rather than state mon- 
itoring, this projection of the set of initial states onto 
a single state resembles a dynamical state preparation 
procedure. The latter enforces a reduction of complex- 
ity in the sense that (after the projection has occurred) 
there is a one-to-one correspondence between the classi- 
cal information contained in the measurement record I(t) 
and the in principle complex quantum information con- 



9 



tained in the state \ip c (t)). In the monitored dynamics, 
this reduction of complexity is manifest in the observed 
transition from quantum to classical properties, i.e., the 
BEC evolves approximately in a coherent state. 

Experimentally, BECs in a double-well potential have 
been thoroughly studied [2E\ |3lJ HE US ES] and the ef- 
fective inter-atomic interaction u is under precise control. 
The investigated unsharp measurement thus represents a 
prototype experiment which would allow to observe the 
transition from an informationally incomplete measure- 
ment (at u = 0) to the informationally complete case 
(at weak inter-atomic interactions) by tuning the inter- 
nal system dynamics. Given the estimated laser intensity 
of the probe field to be 0.75mW for N = 100 bosons, 
these experiments should be realizable with state-of-the- 
art technologies. 

As a future direction of this research line, one may 
explore the possibility to estimate the total number of 



atoms in a BEC placed in a double well by measuring 
only the number of atoms in one well. Furthermore, the 
general concept of quantum to classical transitions, as 
induced by the measurement, should be tested for dif- 
ferent systems. Finally, the full knowledge of the BECs 
state in the presence of dynamics may be employed to 
control the system by unitary feedback depending on the 
measurement record. 
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